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ABSTRACT 

We have studied 14 southern infrared dark clouds (IRDCs) using the data taken from the Mil- 
limetre Astronomy Legacy Team 90 GHz (MALT90) survey and the GLIMPSE and MIPSGAL 
mid-infrared survey of the inner Galaxy. The physical and chemical characteristics of the 14 
IRDCs are investigated using N2H+(l-0), HNC(l-O), HCO+(1-0), and HCN(l-O) molecular 
lines. We find that the 14 IRDCs are in different evolutionary stages from the "starless" to the 
sources with an UCHII region. Three IRDCs are detected to have the star forming activities. 
The integrated intensity ratios /hco+/hcn. ^N2H+/hcn. and /hnc/hcn are all about 1.5, 
which is different from the previous measurements, suggesting that the integrated intensity 
ratios may be affected by the cloud environments. The integrated intensities of HNC, HCO+ 
and HCN show a tight correlation for the 14 IRDCs, implying a close link to the chemical evo- 
lution of these three species in the IRDCs. The derived excitation temperature for each IRDC 
is less than 25 K. We estimate the abundances of the four molecules from 10~^^tolO~^, and 
the average abundance ratios Nhnc/Nhcn = 1-47 ± 0.50, Nhnc/Nhco+ = 1-74 ± 0.22, 
Nhcn/Nhco+ = 1.21 ±0.41. 

Key words: astrochemistry: abundances — ISM: IRDCs — ISM: clouds — stars: formation 
— ISM: molecules-ratio lines 



1 INTRODUCTION 

Infrared dark clouds (IRDCs) are the dark extinction regions of 
high contrast against the bright Galactic mid-infrared background, 
discovered by the Midcourse Space Experiment (MSX ) ( Carey etj 
|al.|1 998| lEg an et al._1998.) and infrared Space Observatory (ISO) 
surveys (Pera ult et al"|1996| ). [Simon et al.| f2006) identified 10,931 
IRDCs candidates using MSX 8.3 /xm data. ^Jackson et al.| ( (2008j ) 
showed that these identified IRDCs are located in the fourth quad- 
rant and first quadrant of the Galaxy at a galactocentric distances of 
6 kpc and 5 kpc, respectively. Additionally, P eretto & Fuller ( 2009 J 
identified 1 1,303 IRDCs candidates using Spitzer 8 iim data. 

Previous researchers suggested that IRDCs were the cold 
(T < 25 K) and dense (10^ cm"^) region s, with a scale of 1^10 
pc and a mass of 10^ ~ 10^ Mq ( |Egan et al.|l 998', 'Care y et al. 
p^98l [2OOOI [Rathborne et al.||2006| ). [Chambers et al.^(200gl ) pro- 



drocarbon (PAH) emission at 8 fim, it is considered to be in the 
finally red core stage. Furthermore, these detected cores have the 
strong dust emission from millimetre and submillimetre bands (Lis 
& Carlstromp994| [Carey et al.. 2000, , Redman et al.,2003, Rath 



posed that the cores within the IRDCs may be in different phase, 
from a quiescent to an active, and finally into a red core. The quies- 
cent cores represent the earliest preprotostellar (starless) core phase 
without infrared signatures, while the active cores have the ex- 
tended and enhanced 4.5 /im emission with an embedded 24 /im 
emission source. When a core shows the polycyclic aromatic hy- 



borne et al.[2005[[Beuther et al.[2005[[Rathborne et al.]2005| [20061 ) 

and only some cores have embedded protostars, indicating that the 
IRDCs may represent the earliest observable stage of high-mass 
star formation. Thus, IRDCs can provide us with an opportunity to 
study the physical and chemical conditions of massive star-forming 
processes in the earliest stage. 

In this paper, we have analyzed 14 southern IRDCs us- 
ing N2H+(l-0), HNC(l-O), HCO+(1-0), and HCN(l-O) molec- 
ular lines from th e Millimetre Astron omy Legacy Team 90 GHz 
(MALT90) survey ( [Foster et al.|20 lit . N2H+, HNC, HCO+, and 
HCN are good tracers of dense gases. And N2H+ is known to be a 
good tracer of the compact center of the cores. Observations of the 
molecular lines could provide valuable information on the physical 
and chemical significance of IRDCs. 



2 DATA 

2.1 IRDCs selection 
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10,931 IRDC candidates. Combined this catalog with the MALT90 
survey, we select 18 southern IRDC candidates, but only 14 of them 
have N2H+(l-0), HNC(l-O), HCO+(1-0), and HCN(l-O) emis- 
sion with high signal-noise ratio. These lines are all good tracers 
of dense gases, but provide slightly different information. N2H^ 
was more resistant to freeze-out on grains than the carbon-bearing 
spec ies (Bergin et al. 2001 ) . HNC was particularly prevalent in cold 
gas ( [Hirota et al.|l998|).HCO^ ofte n showed infall signatures and 
outflow wings ([Rawli ngs et aL|2004[|Fuller et aL|2005| ). The Spitzer 
IRAC 8 /xm data of these selected IRDCs are presented (gray) in 
Figures 1-14, and the basic information are summarized in Tables 
1-2. 



2.2 IRDCs observations 

The data extracted from the Millimetre Astronomy Legacy Team 90 
GHz (MALT90), GLIMPSE and MIPSGAL surveys are analyzed 
toward the 14 southern IRDCs. 

MALT90 is a large international project aimed at character- 
izing high-mass dense cores in the southern sky at 90 GHZ with 
the Mopra 22-m Telescope. The angular resolution of Mopra Tele- 
scope is about 38". The correction for the line intensities to the 
main beam brightness temperature scale is made by using the for- 
mula Tmb = T\/r]u, where r]u is the frequency-dependent beam 
efficiency. The main beam efficiency at 86 GHZ is ?786Ghz = 0.49, 
and a t 110 GHZ is rynocHz 0.44 (| Lo et al.|[255^|Ladd et a l | 
120051 MALT90 data cubes are downloaded from online archiv^ 
The data are reduced by the software CLASS (Continuum and Line 
Analysis Single-Disk Software) and GREG (Grenoble Graphic). 

GLIMPSE is a mid-infrared survey of the inner Galaxy per- 
formed with the Spitzer Space Telescope. MIPSGAL is a survey of 
the same region as GLIMPSE, using the MIPS instrument (24 /xm 
and 70 /am) on Spitzer. We use the mosaicked images of GLIMPSE 
at Spitzer IRAC 8 /xm and MIPSGAL at 24 /xm. S pitzer IRAC 8 fim 
has an angular resolution between 1.5" and 1.9" (|Fazio et al. 2004 
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[Werner et akp004j , and the angular resolution of MIPSGAL 24 fim 
is - 6". 



3 RESULTS 
3.1 Spectra 

Figures 1-14 show the average spectra of N2H+(l-0), HNC(l-O), 
HCO+(1-0), andHCN(l-O) of the 14 southern IRDCs, respectively. 
From each Figure, we see that the N2H^ and HCN lines present 
the hyperfine structure (HFS), and their velocity components blend 
with each other. However, the main velocity component (I23 — O12) 
of N2H+ is detected clearly for most of the IRDCs. The N2H+ and 
HCN lines are fitted using a HFS fit procedure. The fitting results 
are presented in Table 1. From Table 1, we can find that the op- 
tical depths of the N2H^ lines are less than 1 for all the IRDCs, 
indicating that the N2H^ line is optically thin in the IRDCs, which 
agrees with the previous researches. For the HCN line, it is also 
optically thin, which is inconsistent with the previous studies. Con- 
sidering the quality of the HCN lines' data and the unsatisfied HFS 
fits of the HCN lines, it is probable that we underestimate the opti- 
cal depths of the HCN lines. The velocity widths of the N2H^ lines 
are between 1 and 3.2. Calculating the velocity dispersion Ay of 
the optically thin N2H+ line causing by thermal motions: 



http://atoa.atnf.csiro.au/MALT90/ 



where T is the gas kinematic temperature, mobs is the mass of the 
observed species (29 per amu for N2H^), and < m > is the mean 
molecular mass (2.3 per amu). For gas with T = 20 K, all the 
sources have N2H^ line widths AV > AVther ~ 0.68, i.e., hav- 
ing greater nonthermal rather than thermal motions. We assume that 
the nonthermal motions in the N2H^ line width may arise from tur- 
bulence ( [Mardones et al.|1997t . 

In Figures 1-14, the HNC and HCO^ spectra show a wide va- 
riety of line profiles including the double peak, a peak and a "shoul- 
der", a peak skewed to the blue side and single symmetric lines. The 
HNC and HCO^ line shapes differ from source to source but are 
usually similar to each other. The N2H^ line, on the other hand, 
is Gauss toward almost all the sources. In sources with symmetric 
HNC and HCO^ lines, their peak velocity lies very close to that 
of the N2H^ line. In sources with double-peaked HNC and HCO^ 
lines, the N2H^ peak velocity lies between the two peaks (or be- 
tween the peak and the shoulder), indicating that the complex HNC 
and HCO^ line profiles arise from self-absorption at low velocities. 
The Gauss fit is used for the HNC and HCO^ lines and the fitting 
results are listed in Table 2. 

N2H+(l-0), HNC(l-O), HCO+(1-0), and HCN(l-O) all have 
much higher critical density (^ 10^ cm~^) for collisional ex- 
citation. Therefore, there are no other N2H+(l-0), HNC(l-O), 
HCO+(1-0), and HCN(l-O) sources in the line of sight direction. 
For this reason, we can use the Vlsr of the molecules to de- 
termine the kinematic distance to each IRDC, according to the 
rotation curve of |Reid et aL] ( |2009| ), where the Galactic center 
is Ro = 8.4 ± 0.6 kpc and a circular rotation speed is Qo = 
254 ± 16 km s~^. Therefore, four far and four near kinematic dis- 
tances are obtained for each IRDC. The calculating results are in 
Tables 1-2. Since IRDCs are perceived as dark extinction features 
against the Galactic background, it is reasonable to assume that 
all IRDCs are located at the near kinematic distances. Under this 
assumption, we can obtain the average kinematic distance and its 
corresponding error from the four molecules for every IRDC. Here, 
we do not consider the errors resulting from the uncertainties of po- 
sitions and the rotation curve. The final results are listed in Table 
3. 



3.2 The mm and infrared emission in the IRDCs 

Figures 1-14 show also the diagrams of N2H+, HNC, HCO+ and 
HCN integrated intensity superimposing on the Spizter IRAC 8 iim 
and MIPSGAL 24 iim emission images for every IRDC. The in- 
tegrated intensities are calculated for each line in the same veloc- 
ity range presented in Tables 1-2 for each IRDC. Comparing the 
IR emission of each IRDC, 14 IRDCs may be divided into dif- 
ferent evolutional stages, from the "starless " to the sources with 
strong 8 iim emission. From Figures 1-14, we find that the HNC 
andN2H+ emission both match the silhouettes of IRDCs presented 
by 8 iim extinction. Hence, HNC and N2H^ molecules can be used 
to study the morphology of IRDCs in different stages. 

IRDC G003.399-00.399 - In Figure 1, the emission of N2H+ 
and HNC lines show a similar morphology with a single core, but 
which are different from that of HCO^ and HCN. There are 24 iim 
and 8 /xm emission sources close to the peak of N2H^ emission. 

IRDC G003.436-00.572 - In Figure 2, the integrated inten- 
sity maps of N2H+(l-0), HNC(l-O), HCO+(1-0), and HCN(l-O) 
lines all show a morphology extended from south-east to north- 
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Table 1. The physical parameters for molecular lines N2H+ and HCN. And all the parameters are averaged on the pixels 
with the integrated intensity > 5cr and the intensity > 3cr. 



source 


1 


b 


molecular line 


Tex 


^LSR 


Width 




J mb 


Integrated range 




dnear 




deg 


deg 




K 


(kms~l) 


(kms~l) 




(K-kms~l) 


(km s ~ 1 ) 


kpc 


kpc 


G003.399-00.399 


3.310 


-0.398 


N2H + 


19.81±1.62 


6.36 ±0.18 


2.57±0.36 


0.11±0.12 


5.21 ±0.30 


(2,9) 


14.22 


2.55 




3.310 


-0.398 


HCN 


17.71±1.50 


6.22 ±0.28 


2.09±0.55 


0.13±1.30 


2.41 ±0.30 


(2,9) 


14.26 


2.51 


G003.436-00.572 


3.437 


-0.571 


N2H + 


12.05 ±1.27 


2.67 ±0.12 


1.33±0.28 


0.16±0.20 


2.84 ±0.21 


(0,5) 


15.51 


1.26 




3.437 


-0.571 


HCN 


16.28 ±1.24 


2.36 ±0.34 


2.28±0.61 


0.11 ±0.59 


1.72 ±0.22 


(0,5) 


15.64 


1.13 


GOlO.402-00.202 


10.403 


-0.202 


N2H + 


17.22±1.13 


11.70 ±0.13 


1.71±0.29 


0.12±0.10 


3.84 ±0.20 


(8.5,14.5) 


14.78 


1.74 




10.403 


-0.202 


HCN 


14.23 ±1.24 


11.32 ±0.41 


3.09±0.72 


0.12±0.62 


2.28 ±0.23 


(8.5,14.5) 


14.83 


1.70 


GOlO.990-00.083 


10.990 


-0.082 


N2H + 


12.76±1.04 


29.42 ±0.12 


1.96±0.28 


0.15±0.14 


3.74 ±0.19 


(26,32) 


13.22 


3.27 




10.990 


-0.082 


HCN 


14.74±1.03 


28.90 ±0.32 


2.14±0.63 


0.10±0.68 


1.31 ±0.19 


(26,32) 


13.26 


3.24 


G308. 121-00. 152 


308.122 


-0.336 


N2H + 


14.55±1.17 


-47.17 ±0.17 


1.67±0.36 


0.11±0.15 


2.46 ±0.21 


(-50,-44) 


6.71 


3.63 




308.122 


-0.336 


HCN 


18.91±1.07 


-47.21 ±0.24 


2.21 ±0.50 


0.11 ±0.48 


2.03 ±0.20 


(-50,-44) 


6.70 


3.63 


G317.701+00.110 


317.701 


+0.110 


N2H + 


15.65 ±1.40 


-43.59 ±0.13 


1.96±0.31 


0.15±0.15 


3.53 ±0.21 


(-44,-40) 


9.76 


2.63 




317.701 


+0.110 


HCN 


14.53 ±1.46 


-42.03 ±0.33 


1.89±0.51 


0.13±0.65 


1.34 ±0.24 


(-44,-40) 


9.85 


2.54 


G32 1.756+00.029 


321.756 


+0.030 


N2H + 


15.23±1.26 


-32.22 ±0.12 


1.47±0.29 


0.14±0.16 


3.55 ±0.22 


(-35,-29) 


11.21 


1.95 




321.756 


+0.030 


HCN 


12.95 ±1.25 


-32.50 ±0.25 


2.43 ±0.54 


0.17±0.43 


2.76 ±0.23 


(-35,-29) 


11.20 


1.97 


G331.035-00.418 


331.035 


-0.418 


N2H + 


18.94±1.22 


-64.44 ±0.17 


2.72±0.33 


0.11 ±0.06 


6.08 ±0.28 


(-70,-60) 


11.05 


3.63 




331.035 


-0.418 


HCN 


7.41±1.14 


-65.13 ±0.31 


4.16±0.59 


0.38±0.26 


8.11 ±0.30 


(-70,-60) 


11.02 


3.66 


G33 1.708+00.583 


331.708 


+0.584 


N2H + 


23.32±1.17 


-67.46 ±0.15 


3.18±0.28 


0.11 ±0.06 


10.69±0.27 


(-71,-61) 


11.00 


3.77 




331.708 


+0.584 


HCN 


6.94 ±1.14 


-68.09 ±0.33 


2.98±0.61 


0.31±0.71 


3.60 ±0.26 


(-71,-61) 


10.97 


3.80 


G334. 198-00.202 


334.199 


-0.201 


N2H + 


12.65 ±1.04 


-48.03 ±0.15 


1.84±0.38 


0.13±0.15 


3.02 ±0.20 


(-51,-44) 


12.12 


2.98 




334.199 


-0.201 


HCN 


13.68±1.03 


-49.20 ±0.36 


2.00±0.60 


0.10±0.75 


1.37 ±0.20 


(-51,-44) 


12.07 


3.04 


G337.764-00.338 


337.765 


-0.337 


N2H + 


20.69±1.12 


-41.77 ±0.12 


1.92±0.26 


0.11 ±0.07 


4.89 ±0.21 


(-45,-38) 


12.68 


2.85 




337.765 


-0.337 


HCN 


14.53 ±1.07 


-42.46 ±0.23 


1.98±0.50 


0.13±0.30 


2.06 ±0.21 


(-45,-38) 


12.65 


2.88 


G341.942-00.167 


341.944 


-0.166 


N2H + 


22.44±1.35 


-42.40 ±0.15 


2.77±0.29 


0.11 ±0.07 


8.40 ±0.31 


(-47,-37) 


12.78 


3.17 




341.944 


-0.166 


HCN 


12.36±1.38 


-44.01 ±0.28 


3.37±0.55 


0.22±0.37 


6.88 ±0.33 


(-47,-37) 


12.70 


3.25 


G344.726-00.541 


344.727 


-0.540 


N2H + 


10.82±1.14 


-33.15 ±0.14 


1.29±0.33 


0.15±0.33 


2.04 ±0.18 


(-36,-31.5) 


13.27 


2.92 




344.727 


-0.540 


HCN 


15.79±1.07 


-33.13 ±0.29 


1.74±0.58 


0.10±0.68 


1.08 ±0.18 


(-36,-31.5) 


13.27 


2.92 


G345.556+00.026 


345.557 


-0.055 


N2H + 


13.03±1.18 


-16.33 ±0.11 


1.07±0.27 


0.14±0.38 


1.92 ±0.19 


(-19,-14) 


14.48 


1.78 




345.557 


-0.055 


HCN 


15.40±1.16 


-16.64 ±0.32 


2.08±0.56 


0.11 ±0.52 


1.54 ±0.20 


(-19,-14) 


14.45 


1.81 



Table 2. The physical parameters for molecular lines HNC and HCO+. And all the parameters are averaged on the pixels 
with the integrated intensity > 5cr and the intensity > 3cr. 



source 


1 


b 


molecular lines 


Tmb 




Width 


T 


J T^mhdv 


Integrated range 


dfar 


dnear 




deg 


deg 




K 


(kms~l) 


(kms~l) 




(K- kms~l) 


(kms~l) 


kpc 


kpc 


G003.399-00.399 


3.310 


-0.398 


HNC 


1.37±0.34 


6.13 ±0.20 


2.80±0.53 


0.10±0.03 


4.11 ±0.29 


(2,9) 


14.29 


2.48 




3.310 


-0.398 


HCO+ 


1.33±0.35 


6.34 ±0.18 


1.89±0.45 


0.10±0.03 


3.09 ±0.30 


(2,9) 


14.23 


2.54 


G003.436-00.572 


3.437 


-0.571 


HNC 


1.22±0.28 


2.52 ±0.17 


2.10±0.42 


0.10±0.02 


2.71 ±0.21 


(0,5) 


15.57 


1.20 




3.437 


-0.571 


HCO+ 


1.06±0.28 


1.58 ±0.16 


1.47±0.41 


0.08±0.02 


1.73 ±0.21 


(0,5) 


15.98 


0.79 


GO 10.402-00.202 


10.403 


-0.202 


HNC 


1.04±0.25 


11.81 ±0.21 


3.31±0.51 


0.10±0.03 


3.43 ±0.20 


(8.5,14.5) 


14.77 


1.76 




10.403 


-0.202 


HCO+ 


0.98±0.26 


11.93 ±0.26 


4.06±0.61 


0.09±0.03 


3.90 ±0.21 


(8.5,14.5) 


14.75 


1.77 


GOlO.990-00.083 


10.990 


-0.082 


HNC 


0.96±0.23 


28.87 ±0.22 


3.27±0.55 


0.09±0.02 


3.03 ±0.18 


(26,32) 


13.26 


3.23 




10.990 


-0.082 


HCO+ 


0.84±0.23 


28.08 ±0.20 


2.24±0.53 


0.08±0.02 


1.97 ±0.19 


(26,32) 


13.32 


3.18 


G308.121-00.152 


308.122 


-0.336 


HNC 


1.28±0.25 


-47.16 ±0.16 


2.57±0.40 


0.08±0.02 


3.42 ±0.20 


(-50,-44) 


6.71 


3.63 




308.122 


-0.336 


HCO+ 


1.60±0.26 


-47.49 ±0.13 


2.24±0.32 


0.11 ±0.02 


3.99 ±0.21 


(-50,-44) 


6.67 


3.67 


G317.701+00.110 


317.701 


+0.110 


HNC 


1.43 ±0.29 


-42.92 ±0.18 


2.65±0.48 


0.13±0.03 


3.17 ±0.20 


(-44,-40) 


9.80 


2.59 




317.701 


+0.110 


HCO+ 


1.26±0.30 


-42.08 ±0.20 


2.73±0.53 


0.12±0.03 


3.13 ±0.20 


(-44,-40) 


9.85 


2.55 


G32 1.756+00.029 


321.756 


+0.030 


HNC 


1.48 ±0.29 


-32.27 ±0.15 


2.13±0.34 


0.16±0.04 


3.52 ±0.24 


(-35,-29) 


11.21 


1.96 




321.756 


+0.030 


HCO+ 


1.62±0.30 


-32.52 ±0.16 


2.53±0.39 


0.1 8 ±0.04 


4.39 ±0.24 


(-35,-29) 


11.19 


1.97 


G331.035-00.418 


331.035 


-0.418 


HNC 


1.54±0.26 


-64.35 ±0.19 


4.25±0.43 


0.44±0.14 


6.94 ±0.27 


(-70,-60) 


11.05 


3.62 




331.035 


-0.418 


HCO+ 


2.29±0.28 


-64.32 ±0.14 


4.52±0.33 


0.75±0.26 


11.11±0.29 


(-70,-60) 


11.05 


3.62 


G33 1.708+00.583 


331.708 


+0.584 


HNC 


1.23±0.25 


-67.85 ±0.19 


3.23±0.45 


0.38±0.15 


4.37 ±0.26 


(-71,-61) 


10.98 


3.79 




331.708 


+0.584 


HCO+ 


1.33±0.26 


-68.20 ±0.18 


3.02±0.44 


0.42±0.16 


4.44 ±0.27 


(-71,-61) 


10.97 


3.80 


G334. 198-00.202 


334.199 


-0.201 


HNC 


0.92±0.23 


-48.16 ±0.21 


2.87±0.54 


0.09±0.02 


2.82 ±0.20 


(-51,-44) 


12.12 


2.99 




334.199 


-0.201 


HCO+ 


0.88±0.24 


-49.41 ±0.17 


1.58±0.39 


0.09±0.03 


1.89 ±0.21 


(-51,-44) 


12.06 


3.05 


G337.764-00.338 


337.765 


-0.337 


HNC 


1.81±0.26 


-41.97 ±0.12 


2.68±0.30 


0.17±0.03 


5.21 ±0.23 


(-45,-38) 


12.67 


2.86 




337.765 


-0.337 


HCO+ 


1.39±0.26 


-42.72 ±0.16 


2.42±0.38 


0.13±0.03 


3.42 ±0.23 


(-45,-38) 


12.64 


2.90 


G34 1.942-00. 167 


341.944 


-0.166 


HNC 


1.92±0.30 


-43.02 ±0.17 


3.99±0.43 


0.23±0.04 


8.46 ±0.31 


(-47,-37) 


12.75 


3.20 




341.944 


-0.166 


HCO+ 


2.30±0.32 


-43.76 ±0.16 


3.93±0.42 


0.29±0.05 


10.00±0.33 


(-47,-37) 


12.72 


3.24 


G344.726-00.541 


344.727 


-0.540 


HNC 


0.99±0.25 


-32.95 ±0.18 


2.08±0.40 


0.08±0.02 


2.11 ±0.17 


(-36,-31.5) 


13.28 


2.91 




344.727 


-0.540 


HCO+ 


0.94±0.27 


-33.32 ±0.22 


2.47±0.55 


0.08±0.02 


2.18 ±0.19 


(-36,-31.5) 


13.26 


2.93 


G345.556+00.026 


345.557 


-0.055 


HNC 


0.96±0.25 


-16.61 ±0.18 


1.99±0.42 


0.08±0.02 


2.01 ±0.18 


(-19,-14) 


14.45 


1.80 




345.557 


-0.055 


HCO+ 


0.99±0.26 


-16.78 ±0.19 


2.16±0.43 


0.08±0.02 


2.17 ±0.19 


(-19,-14) 


14.44 


1.82 



west. Two compact cores are clearly shown in HCN line emission, 
but we cannot see the obvious IR emission. This IRDC seems to be 
a "starless". 

IRDC GlO.402-00.202 and IRDC GlO.990-00.083 - From 
Figure 3 and Figure 4, we can see that both IRDCs present elon- 
gated structures in all the molecular emission, but have different 
extended directions. For IRDC GlO.402-00.202, the integrated in- 
tensity map of N2H+(l-0) show two cores, while IRDC G10.990- 
00.083 contains three compact cores. At the same time, the spec- 
tral profiles of both IRDCs have double peaks in HNC and HCO^ 
lines, while the optical thin line N2H^(l-0) have a single peak. 
An 24 iim emission source is close to the peak of the molecular 
emission in both IRDCs. 



IRDC G308.121-00.337 - In the diagrams of Figure 5, a north- 
east elongated and compact structure is showed in the four molec- 
ular emission, which contains two cores. An 8 /xm emission source 
is close to the peak of N2H+ emission. The spectra of HNC and 
HCO^ exhibit the asymmetric profile. 

IRDC G317.701-h00.110 - A compact core and an extended 
core are showed in the integrated intensity map of N2H^ and HNC 
lines in Figure 6. At the center of the compact core, there are obvi- 
ous emission at 24 /xm and 8 /xm, here we do not detect the emis- 
sion of the HCN line in this IRDC. 

IRDC G321. 756-^00.029 - Figure 7 shows a compact core 
elongating from SE to NW in the four molecular emission. It 
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Table 3. The integrated intensity ratios of the four species. 



Source 


'^N2H+ /HCN 


^HNC /HCN 


^HCO+ /HCN 


^HCN/HCO+ 


d(kpc) 


G003.399-00.399 


2.16±0.30 


1.71±0.24 


1.28±0.20 


0.78±0.12 


2.52±0.04 


G003 .436-00.572 


1.65±0.24 


1.58±0.24 


1.01±0.18 


0.99±0.18 


1.10±0.30 


GOlO.402-00.202 


1.68±0.19 


1.50±0.18 


1.71±0.20 


0.58±0.07 


1.74±0.04 


GOlO.990-00.083 


2.85±0.44 


2.31±0.36 


1.50±0.26 


0.66±0.12 


3.23±0.05 


G308.121-00.152 


1.21±0.16 


1.68±0.19 


1.97±0.22 


0.51±0.06 


3.64±0.03 


G317.701+00.110 


2.63±0.50 


2.37±0.45 


2.34±0.44 


0.43±0.08 


2.58±0.05 


G321. 756+00.029 


1.29±0.13 


1.28±0.14 


1.59±0.16 


0.63±0.06 


1.96±0.01 


G331.035-00.418 


0.75±0.04 


0.86±0.05 


1.37±0.06 


0.73±0.03 


3.63±0.03 


G331.708+00.583 


2.97±0.23 


1.21±0.11 


1.23±0.12 


0.81±0.08 


3.79±0.02 


G334. 198-00.202 


2.20±0.35 


2.06±0.33 


1.38±0.25 


0.72±0.13 


3.02±0.04 


G337.764-00.338 


2.37±0.26 


2.53±0.28 


1.66±0.20 


0.60±0.07 


2.87±0.03 


G341. 942-00. 167 


1.22±0.07 


1.23±0.07 


1.45±0.08 


0.69±0.04 


3.22±0.05 


G344.726-00.541 


1.89±0.36 


1.95±0.36 


2.02±0.38 


0.50±0.09 


2.92±0.01 


G345.556+00.026 


1.25±0.20 


1.31±0.21 


1.41±0.22 


0.71±0.11 


1.80±0.02 



seems that there are no 24 /xm and 8 /xm emission in the IRDC 
G321. 756-^00.029. 

IRDC G33 1.035-00.41 8 - In Figure 8, the morphologies of the 
four molecular emission all extend from NE to SW. There are only 
one core and an extend structure showed in the N2H^ emission, 
but several cores are identified in other molecular emission. From 
IR emission, we find that it is probable a starless source associ- 
ated with an UCHII region in its southeast ( Bronfman et al.|1996| . 
According to the discussion of section 4.1, it may be an infall can- 
didate. 

IRDC G33 1.708-^00.583 - The four molecular emission 
clearly show two cores in this IRDC in Figure 9. Cyganowski et 
[ar|f200 8 ) identified it as an extended Green Objects (EGO) source 
and an outflow candidate. Here we detect a double-peaked profile 
in HNC and HCO^ lines, implying that it is likely to be an infall 
candidate. 

IRDC G334.198-00.202 - From the diagrams of Figure 10, 
N2H^ emission shows an extended structure associated with the 
24 /j.m and 8 /am emission, while the showed structures are com- 
plicated for other molecular emission because of the weak signal- 
to-noise ratio. 

IRDC G337.764-00.338 - In Figure 11, the four molecular 
emission all show a morphology extending from NE to SW. The 
integrated intensity map of N2H^ line displays a compact core ex- 
tending toward northeast, while the HNC line emission shows that 
IRDC G337.764-00.338 has two cores. The optically thick HNC 
and HCO^ lines present asymmetric profiles. 

IRDC G34 1.942-00. 167 - In Figure 12, the four molecular 
emission show a similar morphology with a single core. |Bronf-| 
|man et al.|(T996l ) suggested that IRDC G341. 942-00. 167 is asso- 
ciated with an UCHII region. We also detect strong 24 /xm and 
8 /xm emission in this cloud. The HNC line profile shows a peak 
and a shoulder, but that of HCO^ is double-peaked, which may be 
caused by an infall motion. 

IRDC G344.726-00.541 and IRDC G345.556-h00.026 - In Fig- 
ures 13-14, the emission of N2H^, HNC and HCO^ lines present 
similar morphology. Seeing the IR emission of them, both of them 
are starless. For IRDC G345.556-h00.026, the HCO+ lines show a 
absorption dip in the redshift relative to the Vlsr- 



4 DISCUSSION 

4.1 infall and outflow 

For self-absorbed optically thick lines, the classical signature of 
infall is a double-peaked profile with the blue shifted peak being 
stronger, or a line asymmetry with the peak skewed to the blue 
side, while optically thin lines should show a single velocity com- 
ponent peaked at the line center. In section 3.2, some IRDCs present 
the blue profiles. In order to provide a strong evidence on whether 
these sources have the infall motion, we plot the map grids of 
HCO+(1-0) for each IRDCs and find only three IRDCs {IRDC 
G331.035-00.418, IRDC G33 1.708-^00.583 and IRDC G341.942- 
00.167) may have the infall motions. Figures 15-16 show the map 
grids towards the three IRDCs, which seems to show the infall fea- 
tures in the whole mapping observations. However, in order to af- 
firm that the spectra in the mappings really show the infall signa- 
tures, we extract the molecular lines from two positions (In Figures 
15-16). In every diagram of the spectra, the optically thin N2H^ 
line is plotted in black color, while the optically thick HNC and 
HCO^ lines are presented in green and blue, respectively. The 
black dash lines and the red dash lines mark the positions of the 
Vlsr of N2H^ line and the absorption dip of optically thick lines, 
respectively. From the mappings and the spectra in Figures 15-16, 
these three IRDCs may have the infall motions in the large-scaled 
regions. We estimate the extent of the infall signature to be up to 
2 X 2armin^ (at least 2.11 pc for IRDC G331.035-00.418, 1.10 
pc for IRDC G33 1.708-^00.583 and 1.40 pc for I RDC G34 1.942- 
00. 167). I n addition, for IRDC G331.708-hOO.583, "Cyganowski et 
[aT] (|2008| ) detected it as an outflow candidate and Yu & Wang 
(|2012) found that it has two outflows, corresponding to the two 
cores in the cloud. At the same time, for IRDC G331.035-00.418 
and IRDC G341. 942-00. 167, we plot their P-V diagrams and do 
not find the outflow signatures. Infall and outflow both are the sign 
of star forming. Hence, we suggest that IRDC G33 1.035-00.41 8, 
IRDC G341. 942-00. 167 and IRDC G33 1.708-^00.583 are forming 
stars, while other IRDCs without the star-forming activity may be 
in much earlier stage. 

4.2 The integrated intensity 

Ratios of the average integrated intensity of N2H^, HNC and 
HCO+ to HCN for each IRDC are presented in Table 3. Consid- 
ering the accuracy of the ratios, we chose the pixels with the inte- 
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grated intensity > 5cr and the intensity > 3a for N2H^, HNC, 
HCO^, and HCN. From Table 3, we find that all /hco+/hcn» 



^N2H+ /HCN 



and Ii 



HNC/HCN 



are almost a constant with around 1.5 



in the error scales for all the IRDCs, implying that the integrated 
intensity ratios seem not to change with the evolution of IRDCs. 
Hsieh et al.| ( |2012| shows that /hcn/hco+ is 1-67 ± 0.83 in the 
starburst Ring and 2.22 ± 0.50 in the Syfert nucleus. It is notice- 
able to us that our /hcn /hco+ of each IRDC is very different from 
those of the starburst Ring and the Syfert nucleus, indicating that 
the integrated intensity ratios probably depend on the cloud envi- 
ronments. 

Figure 17 shows a relationship of the average integrated inten- 
sity between HNC, HCO+ lines and HCN line for the 14 IRDCs. 
We find a tight correlation between HCO+ and HCN (r = 0.98), 
and a linear fitting relationship: 



^HCO+ = (1-32 ± 0.04) X Ihcn + (0.49 ± 0.10) 



(2) 



A very high correlation coefficient r = 0.90 is also found for 
the average integrated intensity of HNC and HCN lines (Figure 17) 
for the 14 IRDCs. The linear fitting result is 



/hnc - (0.77 ± 0.03) X Ihcn + (1.75 =b 0.10) 



(3) 



The results above indicate that there is a close relationship 
for the three species during the process of their chemical evolu- 
tion in the IRDCs. According to this argument, it will contribute 
to determine the dominated chemistry model in the IRDCs through 
the numerical simulation. From Figure 17, we also find that IRDCs 
G331.035-00.418 and G341. 942-00. 167 have larger integrated in- 
tensity of HNC, HCO^, and HCN, which are associated with an 
UCHII region. It seems that the UV radiation field has an influ- 
ence on the chemistry of HNC, HCO^, and HCN molecules in the 
IRDCs, but we need more data to examine this result. 



4.3 Column density and the relationships between the 
abundance ratios and distances 

4. 3. 1 Derivation of physical parameters 

Under an assumption of local thermodynamic equilibrium (LTE), 
the total column density, N, of the molecule can be derived from 
the following formula ( [Scoville et al.|1986] ). 



N = 10^ X 



3k^ 



-exp( 



hz/J Tex + hzy/6k(J + 1) 



4h7r3/i2^2 ^^2kTex' exp(-hz//kTe 



(4) 



where k is the Boltzmann constant, h is the Planck constant, /i is 
the permanent dipole moment of the molecule and J is the rotational 
quantum number of the lower state. Here the permanent dipole mo- 
ment /i of N2H+, HNC, HCO+ an d HCN molecules are 3.40, 3.05, 
3.90, and 2.985H(Mull er et al.|2 0Ql 2005 ), respectively, v is the 
transition frequency, Tex is the excitation temperature, Tmb is the 
main beam brightness temperature which we get from the Gauss 
fits, and Tu is the optical depth. 

Since N2H^ and HCN lines show their hyperfine structures, 
we can obtain their from the HFS fits. The fitted results are listed 
in Table 1 column 8. The main brightness temperature Tmb is given 
by: 



T„,b = f[J(Tex)-J(Tbg)](l-e- 



http://www.astro.uni-koeln.de/cdms/entries 



(5) 



where f is the filling factor, here we assume f = 1 for all the 
molecules, Tbg is the background temperature, and J(T) is defined 
by 

hv 1 



J(T) 



(6) 



k e^p{hiy/kT) 

According to equation (5) and equation (6), we can derive Tex and 
Tu as below: 



k 4(1 - exp(-T.)) 



+ J(Tbg)]-'))"' (7) 



= -lii{l 



f 



-[J(Tex)-J(Tbg)]-^} 



(8) 



So Tex for N2H^ and HCN can be obtained from equation (7). But 
for HNC and HCO^, we assume that they have the same excitation 
temperature to HCN, because they show a tight correlation in their 
integrated intensity. Then the optical depth of HNC and HCO^ 
can be derived from equation (8). Tex and of the four molecules 
are showed in Table 1 and Table 2. Then we can calculate the col- 
umn densities of N2H+, HNC, HCO+ and HCN using equation 
(4), the results are in Table 4. 

4.3.2 The Column density 6)/N2H+, HNC, HCO+, HCN and 
abundance ratios 

Using the equations in section 4.3.1, all the physical parameters 
are derived and presented in Tables 1-4. From Table 1 and Ta- 
ble 2, we find that Tex < 25 K and the optical depths of the 
four molecular lines are less than 1 for all the IRDCs. Consider- 
ing the fact that HNC, HCO+ and HNC fines should be optically 
thick, their optical depths in our results are probably underesti- 
mated. The derived column densities for the four species and the 
abundance ratios Nhnc/Nhcn, Nhco+/Nhcn, Nhnc/Nhco+ 
are listed in Table 4. The uncertainties of the column densities in- 
clude all the errors caused by Tex, opacities and the integrated in- 
tensities. From Table 4, we find that the column densities of the 
four molecules spread over the range of 10^^ ^ 10^^. Then we 
estimate that the abundances of the four molecules using the typ- 
ical H2 column density (0.9 ~ 4.6) x 10^^ cm~^ in southern 
infrared dark clouds ( [Vasyunina et al.|[2009| ) are in the range of 
10~^^ ^ 10~^ , wh ich are consistent with the results of Vasyun- 
ina etal.V201l') and'Zinchenko et al.'(2009) for N2H+, HNC and 



HCN molecules, but ~ one order lower for HCO^ molecule, even 
over two order lower than the results of |Sanhueza et al.| ( |2012| ). 
This abundant difference of HCO^ molecule is mainly caused by 
the optical depth. HCO^ line is optically thin in our study. From 
the abundance ratios of HNC, HCO+, and HCN for each IRDC 
in Table 4, we find that HNC is more abundant than HCO^ and 
HCN, except for IRDC G33 1.035-00.41 8, which agree with the 
fact that HNC molecule is the tracer of the cold gas. And the de- 
rived Nhnc/Nhcn = 1.47 ± 0.50 is consistent with the results 



in the dark cloud cores (Hirota et al.j 1998 ). Nhcn /Nhco+ is ap- 



proximately equal to 1, implying the similar origin and chemistry 
evolution of these two molecules. Furthermore, we also calculate 
the average abundance ratios (Nhnc/Nhcn, Nhnc/Nhco+ 
Nhcn /Nhco+) of the 14 southern IRDCs, which are presented in 
Table 5. From Table 5, we find that Nhcn/Nhco+ is almost the 
same in the error scale in the three different environments, sug- 
gesting that the abundance ratio of HCN to HCO^ may be not 
affected by the environments. The differences of other two ratios 
Nhnc /Nhcn and Nhnc /Nhco+ indicate that IRDCs may repre- 
sent the chemistry of earlier star formation. 
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Table 4. The column density of N2H+, HNC, HCO+ and HCN and the abundance ratios of HNC, HCO+ and HCN. 



Source 


%H + 


Nhnc 


^HCO + 


Nhcn 


Nrnc/Nhcn 


Nhco+/^hcn 


Nhnc/Nhco+ 




(Xl0l2)cm-2 




(Xl0l2)cm-2 










G003. 399-00. 399 


11.40±1.15 


9.27±0.87 


4.38±0.51 


6.29 ±4.03 


1.47 ±0.95 


0.70±0.45 


2.12±0.32 


G003.436-00.572 


4.90 ±0.64 


5.76±0.54 


2.29±0.30 


4.16 ±1.32 


1.38 ±0.46 


0.55±0.19 


2.52±0.41 


GOlO.402-00.202 


7.75 ±0.65 


6.66±0.55 


4.73±0.38 


5.08 ±1.63 


1.31±0.43 


0.93±0.31 


1.41±0.16 


GOlO.990-00.083 


6.55 ±0.62 


5.99±0.45 


2.44±0.26 


2.93 ±1.07 


2.04±0.76 


0.83±0.32 


2.45 ±0.32 


G308. 121-00. 152 


4.35 ±0.53 


8.01±0.58 


5.96±0.41 


5.47 ±1.40 


1.46±0.39 


1.09±0.29 


1.34±0.13 


G317.701-h00.110 


7.03 ±0.76 


6.34±0.58 


3.91±0.36 


3.05 ±1.12 


2.08±0.79 


1.28 ±0.49 


1.62±0.21 


G321.756h-00.029 


6.83 ±0.75 


6.64±0.61 


5.24±0.44 


6.04 ±1.38 


1.10±0.27 


0.87±0.21 


1.27±0.16 


G331.035-00.418 


12.90±0.92 


11.40±1.03 


13.10±1.68 


16.10±2.12 


0.71±0.11 


0.81±0.15 


0.87±0.14 


G331.708-hOO.583 


26.50±1.44 


6.87±0.71 


4.43±0.48 


6.60 ±2.20 


1.04±0.36 


0.67±0.24 


1.55±0.23 


G334. 198-00.202 


5.08 ±0.54 


5.32±0.46 


2.24±0.27 


2.92 ±1.15 


1.82±0.73 


0.77±0.32 


2.38±0.35 


G337.764-00.338 


11.10±0.75 


10.60±0.71 


4.29±0.36 


4.69 ±0.86 


2.26±0.44 


0.91±0.18 


2.47 ±0.27 


G341. 942-00. 167 


20.20±1.36 


16.00±1.25 


12.20±0.96 


15.30±2.93 


1.05 ±0.22 


0.80±0.17 


1.31±0.15 


G344.726-00.541 


3.26 ±0.59 


4.35±0.41 


2.82±0.28 


2.53 ±0.94 


1.72±0.66 


1.11 ±0.43 


1.54±0.21 


G345.556-h00.026 


3.35 ±0.69 


4.07±0.42 


2.76±0.28 


3.58 ±1.03 


1.14±0.35 


0.77±0.24 


1.47±0.21 



Table 5. The average abundance ratios among different environments 



Environment 


Nhnc/Nhcn 


Nhnc/Nhco+ 


Nhcn/Nhco+ 


References 


IRDCS 
ISM 

Star forming regions 


1.47±0.50 
0.21±0.05 
0.21±0.06 


1.74±0.23 
0.5±0.3 


1.21±0.41 
1.47±0.86 
1.9±0.9 


Present work 
'Liszt & Lucas'j200l| 
Godard et al. (20f^ 



4.3.3 The relationship between abundance ratios and distances 

Figure 18 shows the relationships between the average abundance 
ratios and the distances of the 14 southern IRDCs. Nhco+ /Nhcn, 
Nhco+ /Nhnc and Nhcn/Nhnc all show a linear relationship 
with the distances. The relationships are: 

Nhco+/Nhcn = (0.06 ± 0.07) x d + (0.66 db 0.21) (9) 

Nhco+/Nhnc = (0.03 ± 0.03) x d + (0.47 ± 0.07) (10) 

Nhcn /Nhnc = (0.03 ± 0.08) x d + (0.56 =b 0.23) (1 1) 

From above relationships, we find the abundance ratios 0/N in- 
creasing slowly with the distances of the IRDCs, and a small in- 
crease is also found for the abundance ratio Nhcn/Nhnc to the 
distance. Since the errors and limit data, the more studies should be 
done to check this conclusion. 



5 SUMMARY 

We do the research of 14 southern IRDCs with N2H+(l-0), 
HNC(l-O), HCO+(1-0), and HCN(l-O) lines of the MALT90 survey 
and Spitze 8 /xm, 24 /xm data. The integrated intensity diagrams of 
the four molecular lines are mapped and the physical parameters 
are obtained for the 14 southern IRDCs. We also discuss the kinetic 
processes and explore the chemical features. Our main results are 
summarized as follows. 

1. The Spitzer images show that the 14 IRDCs are in different 
evolutional stages, from "starless cores" with no IR emissions to 
"red cores" with the strong 8 /xm emission. 

2. According to the optically thick HNC and HCO+ line 
profiles and mappings, three IRDCs (IRDC G331.035-00.418, 
G33 1.708-^00.583, and G341. 942-00. 167) are found to have the in- 
fall motions, while other IRDCs without the star-forming activity 
may be in much earlier stage. 

3. The integrated intensity of HNC, HCO+ and HCN correlate 
well with each other for the 14 IRDCs, implying a close link to their 
chemistry evolution in the IRDCs. 

4. The obtained physical parameters show that Tex is < 25 K 



for all the 14 IRDCs and N2H+ line is optically thin in the IRDCs. 
The column densities of the four molecules span up to two orders 
10^^ ^ 10^^ and their corresponding abundances are in the range 
of 10~^^ to 10~^. The average abundance ratios Nhnc /Nhcn = 
1.47 =b 0.50, Nhnc /Nhco+ = 1-74 ± 0.22 show significant dif- 
ference with those of ISM and star forming regions, indicating that 
the environment of earlier star formation may be different. How- 
ever, Nhcn/Nhco+ — 1-24 ± 0.41 is almost the same with that 
of the other environments, suggesting that the abundance ratio of 
HCN to HCO^ may be not affected by the environment. 
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Figure 1. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC G003.309-00.399. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /xm. The red dash line represents the Vlsr of N2H+ line 
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Figure 2. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC G003.436-00.572. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash line represents the Vlsr of N2H+ line 
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Figure 3. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC GOlO.402-00.202. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash line represents the Vlsr of N2H+ line 
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Figure 4. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC GOlO.990-00.083. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash line represents the Vlsr of N2H+ line 
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Figure 5. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC G308. 121-00.152. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash line represents the Vlsr of N2H+ line 
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Figure 6. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC G3 17.701+00.1 10. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash line represents the Vlsr of N2H+ line 
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Figure 7. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC G321. 756+00.029. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash line represents the Vlsr of N2H+ line 
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Figure 8. The spectra and the integrated intensity maps of N2 H+ , HNC, HCO+ , and HCN in IRDC G33 1 .035-00.4 1 8. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash line represents the Vlsr of N2H+ line 
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Figure 9. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC G331. 708+00.583. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash line represents the Vlsr of N2H+ line 
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Figure 10. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC G334. 198-00.202. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash line represents the Vlsr of N2H+ line 
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Figure 11. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC G337.764+00.338. The spectra are in the left; the 
grayscal in the middle plane is Spitzer IRAC 8 /j,m and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash Hne represents the Vlsr of N2H+ 
line 
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Figure 12. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC G341.942-00.167. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash line represents the Vlsr of N2H+ line 
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Figure 13. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC G344.726-00.541. The spectra are in the left; the grayscal 
in the middle plane is Spitzer IRAC 8 /im and that in the right plane is Spitzer MIPSGAL 24 /im. The red dash line represents the Vlsr of N2H+ line 
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Figure 14. The spectra and the integrated intensity maps of N2H+, HNC, HCO+, and HCN in IRDC G345.556+00.026. The spectra are in the left; the 
grayscal in the middle plane is Spitzer IRAC 8 /j,m and that in the right plane is Spitzer MIPSGAL 24 /j,m. The red dash Hne represents the Vlsr of N2H+ 
line 
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Figure 15. The upper plane: the HCO+ map grid of IRDC G331.035-00.418 and the extracted spectra of N2H+ (black), HNC(green) and HCO+(blue) lines 
in two positions, corresponding to the red lines in the map grid. The black dash line and the red dash line mark the position of the Vlsr of N2H+ line and 
the absorption dip of the optically thick lines; The bottom plan: the HCO+ map grid of IRDC G341. 942-00. 167 and the extracted spectra of N2H+(black), 
HNC(green) and HCO+(blue) lines in two positions, corresponding to the red lines in the map grid. The black dash line and the red dash line mark the position 
of the Vlsr of N2H+ line and the absorption dip of the optically thick lines. 



© 2002 RAS, MNRAS 000,[T]-?? 



A Systemic Study of 14 Southern Infrared Dark Clouds with the N2H+, HNC, HCO^, and HCN Lines 23 



IRDC G331 .708+00.583 



0,5 



-0.5 



. w^^A^H i*v>yW' "STtAtv* 



(-0.3, -0.15) 



NjH"" (1-0) 
HNC (1-0) 

HCO'' (1-0) 



-90 -80 -70 -60 

Velocity (km/s) 



(0.3,-0.15) 

id 


! NjH"" (1-0) 
HCO"" (1-0) 







0.5 



k farcminl 



-0.5 



-90 -80 -70 -60 

Velocity (km/s) 



Figure 16. the HCO+ map grid of IRDC G33 1.708+00.583 and the extracted spectra of N2H+ (black), HNC(green) and HCO+(blue) lines in two positions, 
corresponding to the red lines in the map grid. The black dash line and the red dash line mark the position of the Vlsr of N2H+ line and the absorption dip 
of the optically thick Hnes. 
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Figure 17. Left: plot averaged integrated intensity of HCO+ vs. those of HCN for the 14 IRDCs. Right: same for HNC vs. HCN. The red dash lines represent 
the linear fitting results. 
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Figure 18. Left: plot average abundance ratios of HCO+ to HCN vs. the distances of the 14 IRDCs. Middle: plot average abundance ratios of HCO+ to HNC 
vs. the distances of the 14 IRDCs. Right: plot average abundance ratios of HCN to HNC vs. the distances of the 14 IRDCs. The red dash lines represent the 
linear fitting results. 
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